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ABSTRACT 

We study the effects of very strong magnetic fields on the equation of 
state (EOS) in multicomponent, interacting matter by developing a covariant 
description for the inclusion of the anomalous magnetic moments of nucleons. 
For the description of neutron star matter, we employ a field-theoretical 
approach which permits the study of several models which differ in their 
behavior at high density. Effects of Landau quantization in ultra-strong 
magnetic fields {B > 10^^ Gauss) lead to a reduction in the electron chemical 
potential and a substantial increase in the proton fraction. We find the generic 
result for B > 10^* Gauss that the softening of the EOS caused by Landau 
quantization is overwhelmed by stiffening due to the incorporation of the 
anomalous magnetic moments of the nucleons. In addition, the neutrons become 
completely spin polarized. The inclusion of ultra-strong magnetic fields leads 
to a dramatic increase in the proton fraction, with consequences for the direct 
Urea process and neutron star cooling. The magnetization of the matter never 
appears to become very large, as the value of \H/B\ never deviates from unity 
by more than a few percent. Our findings have implications for the structure of 
neutron stars in the presence of large frozen-in magnetic fields. 
Subject headings: stars: neutron stars - equation of state - stars: magnetic fields 
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1. INTRODUCTION 

Recent observational and theoretical studies motivate the investigation of the effects 
of ultra-strong magnetic fields {B > 10^^ Gauss) on neutron stars. Several independent 
arguments link the class of soft 7— ray repeaters and perhaps certain anomalous X-ray 
pulsars with neutron stars having ultra strong magnetic fields - the so-called magnetars 
(Paczyhski 1992; Thompson & Duncan 1995, 1996; Melatos 1999). In addition, two of the 
four known soft 7— ray repeaters directly imply, from their periods and spin-down rates, 
surface fields in the range 2 — 8 x 10^^ Gauss. Kouvehotou et al. (1998, 1999) argue from 
the population statistics of soft 7— ray repeaters that magnetars constitute about 10% of 
the neutron star population. While some observed white dwarfs have large enough fields to 
give ultra-strong neutron star magnetic fields through flux conservaton, it does not appear 
likely that such isolated examples could account for a significant fraction of ultra-strong 
field neutron stars. Therefore, an alternative mechanism seems necessary for the creation 
ultra-strong magnetic fields in neutron stars. Duncan & Thompson (1992, 1996) suggested 
that large fields (up to 3 x 10^^ x (1 ms/Pi) Gauss, where Pj is the initial rotation period) 
can be generated in nascent neutron stars through the smoothing of differential rotation 
and convection. 

These developments raise the intriguing questions: 

(1) What is the largest frozen-in magnetic field a stationary neutron star can sustain?, and, 

(2) What is the effect of such ultra-strong magnetic fields on the maximum neutron star 
mass? 

The answers to both of these questions hinge upon the effects strong magnetic fields 
have both on the equation of state (EOS) of neutron-star matter and on the structure of 
neutron stars. In this paper, we will focus on the effects of strong magnetic fields on the 
EOS. Subsequent work will be devoted to investigating the effects of strong fields on the 
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structure of neutron stars, incorporating the EOSs developed in this work. 

The magnitude of the magnetic field strength B needed to dramatically affect neutron 
star structure directly can be estimated with a dimensional analysis (Lai & Shapiro 1991) 
equating the magnetic field energy ~ (47ri?^/3)(i?^/87r) with the gravitational binding 
energy Eb.e. ^ GM^/R, yielding 5 - 2 x lO^^ (M/1.4Mo) (i?/10 km)"^ Gauss, where M 
and it! are, respectively, the neutron star mass and radius. 

The magnitude of B required to directly influence the EOS can be estimated by 
considering its effects on charged particles. Charge neutral, beta-equilibrated, neutron-star 
matter contains both negatively charged leptons (electrons and muons) and positively 
charged protons. Magnetic fields quantize the orbital motion (Landau quantization) of these 
charged particles. Relativistic effects become important when the particle's cyclotron energy 
eh, B/{mc) is comparable to it's mass (times c^). The magnitudes of the so-called critical 
fields are B^ = (hc/e) Ae"^= 4.414 x 10^=^ Gauss and B^ = {mp/mefB^ = 1.487 x 10^° 
Gauss for the electron and proton, respectively ( Ae= h/rrieC ~ 386 fm is the Compton 
wavelength of the electron) . It will be convenient to measure the field strength B in units of 
Sg, viz., B* = B/B^. When the Fermi energy of the proton becomes significantly affected 
by the magnetic field, the composition of matter in beta equilibrum is significantly affected. 
In turn, the pressure of matter is significantly affected. We show that this occurs when 
B* ~ 10^, and will lead to a general softening of the EOS. 

In neutron stars, magnetic fields may well vary in strength from the core to the surface. 
The scale lengths of such variations are, however, usually much larger than the microscopic 
magnetic scale Im-, which depends on the magnitude of B. For low fields, for which the 
quasi- classical approximation holds, Im — (Ae^ /5*)(37r^ne)^/^ ~ 10^(ne/ns)^/^/5* fm, 
where rig is the number density of electrons and Ug is the normal nuclear saturation 
density (about 0.16 fm~^). For high fields, when only a few Landau levels are occupied. 
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~ 27r^ne(Ae^ /B*y ^ 7 x 10^{ne/ns)/B*^ fm. In either case, the requirement that 
R >> Im is amply satisfied; hence, the magnetic field B may be assumed to be locally 
constant and uniform as far as effects on the EOS are concerned. 

In non-magnetic neutron stars, the pressure of matter ranges from 2 — 5 MeV fm~^ 
at nuclear density to 200 — 600 MeV fm~'^ at the central density of the maximum 
mass configuration, depending on the EOS (Prakash et al. 1997). These values may 
be contrasted with the energy density and pressure from the electromagnetic field: 
Ef = Pf = 5V(87r) = 4.814 X lO'^B* ^ MeV fm"^. The field contributions can dominate 
the matter pressure for > 10^ at nuclear densities and for B* > 10^ at the central 
densities of neutron stars, and must therefore be included whenever the field dramatically 
infiuences the star's composition and matter pressure. 

In strong magnetic fields, contributions from the anomalous magnetic moments of 
the nucleons must also be considered. Experimentally, Hp = ji^ (fi'p/2 — 1) for the proton, 
i^n — I^Ngnf^ for the neutron, where /xjv is the nuclear magneton and gp = 5.58 and 
Qn — —3.82 are the Lande g-factors for the proton and neutron, respectively. The energy 
\K>n + K'plB ~ 1.67 X 10~^B* MeV measures the changes in the beta equihbrium condition 
and to the baryon Fermi energies. Since the Fermi energies range from a few MeV to tens of 
MeV for the densities of interest, it is clear that contributions from the anomalous magnetic 
moments also become significant for B* > 10^. We demonstrate that for such fields, 
complete spin polarization of the neutrons occurs, which results in an overall stiffening of 
the EOS that overwhelms the softening induced by Landau quantization. 

In magnetized matter, the stress energy tensor contains terms proportional to HB, 
where H = B + Ai^Ai and M. is the magnetization (Landau, Lifshitz & Pitaevskii 1984). 
Thus, extra terms, in addition to the usual ones proportional to S^, are introduced into 
the structure equations (Cardall et al. 1999). The magnetization in a single component 
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electron gas has been studied extensively (Blandford & Hernquist 1982) for neutron star 
crust matter. We generalize this formulation to the case of interacting multicomponent 
matter with and without the effects of the anomalous magnetic moments. We find that 
deviations of H from B occur for field strengths B* ^ 10^. 

Although the effects of magnetic fields on the EOS at low densities, relevant for 
neutron star crusts, has been extensively studied (see for example, Canuto & Ventura 1977; 
Pushiki, Gudmundsson & Pethick 1989; Fushiki et al. 1992; Abrahams & Shapiro 1991; 
Lai & Shapiro 1991; Rognvaldsson et al. 1993, Thorlofsson et al. 1998), only a handful 
of previous works have considered the effects of very large magnetic fields on the EOS of 
dense neutron star matter (Chakrabarty 1996; Chakrabarty, Bandyopadhyay, & Pal 1997, 
Yuan & Zhang 1999). Lai and Shapiro (1991) considered non-interacting, charge neutral, 
beta-equilibrated matter at subsaturation densities, while Chakrabarty and co-authors 
studied dense matter including interactions using a field-theoretical description. These 
authors found large compositional changes in matter induced by ultra-strong magnetic fields 
due to the quantization of orbital motion. Acting in concert with the nuclear symmetry 
energy. Landau quantization substantially increases the concentration of protons compared 
to the field-free case, which in turn leads to a softening of the EOS. This lowers the 
maximum mass relative to the field-free value. In these works, however, the electromagnetic 
field energy density and pressure, which tend to stiffen the EOS, were not included. 
In addition, changes in the general relativistic structure induced by the magnetic fields 
(studied in detail by Bocquet et al. 1995 who, however, omitted the compositional changes 
in the EOS due to Landau quantization) were also ignored. Thus, the combined effects 
of the magnetic fields on the EOS and on the general relativistic structure remain to be 
determined. 

Compared to these earlier works, we make several improvements in the calculation of 
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the EOS. These improvements include (1) a study of a larger class of field-theoretical models 
in order to extract the generic trends induced by Landau quantization, (2) the development 
of a covariant description for the inclusion of the anomalous magnetic moments of the 
nucleons, and (3) a detailed study of magnetization of interacting multicomponent matter 
with and without the inclusion of the anomalous magnetic moments. We also provide 
simple analytical estimates of when each of these effects begin to significantly influence the 
EOS. Our future work will employ the EOSs developed in this work to complete a fully 
self-consistent calculation of neutron star structure including the combined effects of the 
direct effects of magnetic fields on the EOS and general relativistic structure. 

In §2, we present the field-theoretical description of dense neutron star matter including 
the effects of Landau quantization and the nucleon anomalous magnetic moments. Section 
3 contains a detailed study of the effects of Landau quantization on the EOS for two classes 
of Lagrangians. In addition to providing contrasts with earlier work, our results highlight 
the extent to which the underlying interactions affect the basic findings. This section 
also includes new theoretical developments concerning the magnetization of interacting, 
multicomponent matter. Section 4 is devoted to the effects of the anomalous magnetic 
moments on the EOS. Here results for a charge neutral neutron, proton, electron, and 
muon gas are compared with those for interacting matter to asses the generic trends. Our 
conclusions and outlook, including the possible effects of additional components such as 
hyperons, Bose condensates and quarks, are presented in §5. The covariant description 
for the inclusion of the anomalous magnetic moments of the nucleons is presented in the 
Appendix, where explicit formulae for the nucleon Dirac spinors and energy spectra are 
derived. Except where necessary, we use units wherin h and c are set to unity. 
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2. THEORETICAL FRAMEWORK 

For the description of the EOS of neutron-star matter, we employ a field-theoretical 
approach in which the baryons (neutrons, n, and protons, p) interact via the exchange of 
a — LU — p mesons. We study two classes of models, which differ in their behavior at high 
density. The Lagrangian densities associated with these two classes are (Boguta & Bodmer 
1977, Zimanyi & Moszkowski 1990) 

jCi = A - f 1 - ^) %mb^i, + jCm + jCi, 

\ rrib / 

Cn = (l + ^)Cb-%rrn,^b + Cm + Ci. (1) 
V nib J 

The baryon (6 — n,p), lepton {I — e,ii), and meson {a, cu, and p) Lagrangians are given by 
Cb = + qb^f.A^' - g^.^^.uj" - gp,n,l;,p^ - Hh(T^^i.F^''') *6 , 

Cm = \d,ad^a-]^mla^-U{a) + ]^niluj,uj^-\n^^n,, 

+ \mlp,ff - Ip'^'P.u - -/^"F,. , (2) 

where ^b a-nd ■0/ are the baryon and lepton Dirac fields, respectively. The nucleon 
mass and the isospin projection are denoted by nib and rs^, respectively. The mesonic 
and electromagnetic field strength tensors are given by their usual expressions: 

= d^uj^ - d^uj^, P^^ = d^p^ - d^p^^, and F^^ = d^A^ - d^A^,. The strong interaction 
couplings are denoted by g, the electromagnetic couplings by q, and the meson masses by 
m all with appropriate subscripts. The anomalous magnetic moments are introduced via 
the coupling of the baryons to the electromagnetic field tensor with o"^,^ = 2 [^/"' ^^'\ 
strength Kb- We will contrast results for cases with Kb — and Kb taken to be their measured 
values. The quantity U (a) denotes possible scalar self-interactions. It is straightforward to 
include self interactions between both the vector uj and the iso- vector p mesons (Miiller & 
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Serot 1996). Although the electromagnetic field is included in Cj and £77, it assumed to 
be externally generated (and thus has no associated field equation) and only frozen-field 
configurations will be considered. 

The thermodynamic quantities will be evaluated in the mean field approximation, in 
which the mesonic fields are assumed to be constant. The field equations are 

'^'(^) + ^^ = <l .^..2 (3) 



^b9a, (^) nl for Cn , 



'^'(^°) =E^?-.^'" (4) 



m 



Hp°) = Y.9p,r3^nb, (5) 



a- {p - qiA) + (3mi] i^i = E^i , (6) 
a- (p- qbA) + f3ml\ = {E - - (?p,r3,p°)^b , (7) 



where the effective baryon masses are 



mi 



l-i^) for £7 

nit ' 



and nl is the scalar number density. The scalar self-interaction is taken to be of the form 
(Boguta & Bodmer 1977; Glendenning 1982, 1985) 

U{a) = Um^^{g,,af + \c{g,,aY , (9) 

where the in the first term is included to make h dimensionless. In charge neutral, beta 
equilibrated matter, the conditions 

f^n - f^p = fJ'e = f^fi , (10) 

np = ne + n^, (11) 

also apply. Given the nucleon-meson coupling constants and the coefficients in the scalar 
self-interaction, equations (^ through (|TT]) may be solved self consistently for the chemical 
potentials, /Xj, and the field strengths, a, tu^, and p°. 
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3. EFFECTS OF LANDAU QUANTIZATION 



From equation d), the energy spectrum for the leptons is (see, for example, Canuto & 
Ventura 1977) 

Ei = 



where 



m. 



I 2 



(12) 



(13) 



Here, n is the principal quantum number and cr^ (not to be confused with the scalar field a) 
is the spin along the magnetic field axis, is the component of the momentum along the 
magnetic field axis. The quantity v = n + 1/2 — {l/2){qi/\qi\)(Jz characterizes the so-called 
Landau level. Equation (^ gives the energy spectrum for the protons as 



Ep = \jkl + mix + QLOpUJ^ - 9pp^P° , 



(14) 



where m^^^ is obtained by replacing mi on the right hand side of equation (|l^) by m*. The 
neutron energy spectrum is that of the free Dirac particle, but with shifts arising from the 



scalar, vector, and isovector interactions: 



(15) 



At zero temperature and in the presence of a constant magnetic field B, the number 
and energy densities of charged particles are given by 



ni=i. 



I /I I R ^max 

I y i I \ ^ \ ^ 



47r2 



CTz n=0 . 



E\k\ 



n=0 



Ef + kin,.. 



ml 



(16) 



(17) 



Above, fcj „ is the Fermi momentum for the level with the principal quantum number n 
and spin and is given by 



Ef -^n,a. 
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The summation in equation (16) is terminated at Umax, which is the integer preceeding 
the value of n for which /cj ^ is negative. The Fermi energies are fixed by the chemical 
potentials 



E) = lib- Qu^^uJ 







(19) 

9pj3,p'. (20) 
For the protons, the scalar number density may be determined to be (Chakrabarty 1996) 



\qp\Bml 
27r2 



EE In 



E^f + F, 



(21) 



o"z n=0 

The number, energy, and scalar number densities of the neutrons are unchanged in form 



from the field-free case 



fir. 



j^n 3 

3^ 



27r2 



Efkf 



ml 



2E] ^k] 



mz 



(22) 
(23) 

(24) 



The total energy density of the system is 



111 

-mlul + -mlpl + -mla^ + [/(a) 



(25) 



where the last term is the contribution from the electromagnetic field. Use of equations 
(10) and (11), which are satisfied in charge neutral beta-equilibrated matter, in the general 
expression for the pressure, P = J2i fJ'i^i ~ e {i = n,p, e, and /i), allows the pressure to be 
written only in terms of the neutron chemical potential through the relation P = finTib — £■ 
In fact, utilizing the appropriate relations satisfied by the various chemical potentials 
and the number densities involved in the charge neutrality condition, it is easily verified 
that this relation is satisfied even in the presence of additional components such as 
strangeness-bearing hyperons, Bose condensates (pion or kaon), and quarks, which may 
likely exist in dense neutron-star matter. 
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3.1. Magnetization 



The magnetic field strength, H, is related to the energy density by (Landau, Lifshitz & 
Pitaevskii 1984) 

ii' = 47r(||j =B + 47rM, (26) 

where M. is the magnetization. This is equivalent to the set of equations 

drib _ 
dB ~ ' 

^ dsi dsi dE) dsp dsp dE^f dSn de^ dE] 

dB dE) dB dB dEF^ dB dB dE] dB 

9 ndp^ 9 r^doj^ de da , 

The first of these gives 

drin dEJ dn,n drip drip dE-j 



dE] dB dB dB dE^f dB ' 
Using the conditions of charge neutrality and chemical equilibrium, one has 



From the field equations and the definition of the scalar density, 

-0, 

_ gp drip 





9u 


drib 


dB 


ml 


dB 


dp^ 


9p_ 


1 d 


dB 


ml 


2dB 


de 


= 0. 








d^ 





m^ dB 



Note also that 



dE} f dE) 
Utilizing these results, equation (26) becomes 



(28) 



dE) dB dB dB dE^f dB ' ^ ' 



(30) 



(31) 



' (32) 
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where 



Ti 



dsi 
OB 



dm 



^ OB dB ^ dB 
drin 



+ 



dsp 
dB 



II 



E\ + E}-E]-g,p\ 



(33) 



Note that chemical equihbrium ensures that Tn — whence the magnetization takes the 
general form 

Z i^-B}^- (34) 



i=e,iJ,,p,n \ / 

In the case under current consideration, inserting the explicit forms of the energy density 



and number density yields the result 



Ei - Elm B ^ 1 1 , , 



m. 



n,az 



(35) 



This result generahzes the result of Blandford & Hernquist (1982) for an electron gas to the 
case of a multi-component system including interacting nucleons. That the functional form 
of A4 for interacting nucleons is the same as that for non-interacting particles stems from 
the fact that, in the mean field approximation, the field equations for the nucleons reduces 
to the Dirac equation for a free particle, but with an effective mass m*. 



3.2. Results 

In Table 1, we list the various nucleon- meson and meson self-interaction couplings for 
the two classes of models chosen for this study. In each case, the couplings were chosen 
to reproduce commonly accepted values of the equilibrium nuclear matter properties: the 
binding energy per particle B/A, the saturation density Ug, the Dirac effective mass m'^/mn, 
the compression modulus Kq, and the symmetry energy a^y^. The high-density behavior of 
the EOS is sensitive to the strength of the meson couplings employed and the models chosen 
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encompass a fairly wide range of variation. The HS81 model, which has a rather high 
compression modulus, allows us to contrast our results with those of Chakrabarty (1996) 
who also used HS81 in the case when Landau quantization is considered, and to assess 
the effects of the inclusion of magnetic moments. Models HS81 and GM1-GM3 employ 
linear scalar couplings {Cj), while the ZM model employs a nonlinear scalar coupling (Cu), 
which is reflected in the high density behaviors of m*^/mn- Thus, comparison of the HS81, 
GM1-GM3 and ZM models allows us to contrast the effects of the underlying EOS. 

In Figure |l], we show results of some physical quantitites of interest for our baseline 
case, model GM3. At supernuclear densities and in the absence of a magnetic field, the 
matter pressure is dominated by the baryons principally due to the repulsive nature of the 
strong interactions. Even up to the central density in a neutron star, the proton fraction 
remains sufficiently small that the neutrons dominate the total pressure. 

The magnitude of the magnetic field B required to induce significant changes in 
the EOS may be estimated in a straightforward manner. In the presence of a magnetic 
field, the contributions from the protons become significant when only one Landau level 
is occupied, i.e., when the protons are completely spin polarized. This happens when 
qpB/Tic > (27r^np)^/^. Therefore, we arrive at the estimate B* >X^ (27r''y^^n^)^/^ for 
quantum effects to dominate. The proton fraction Yp = Up/ub, which depends upon both 
the density and the magnetic field, typically lies in the range 0.1-0.7. As a result the term 
in parentheses is of order 1 fm~^ for densities > O.ln^, and Ag^— 1.5 x 10^ fm^. Thus, the 
magnetic field necessary to introduce significant contributions from the protons is of order 
B* ~ 10^, which is well below the proton critical field B^ = {mp/me)'^Bl = 1.49 x 10^° 
Gauss (or i?* = 3 x 10^) for which protons begin to become relativistic. 

The results in Figure |I] were obtained by accounting for all of the allowed Landau 
levels. Indeed, we notice that the matter pressure Pm, the effective mass m*, and the 
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concentrations 1^ = ni/rib begin to differ significantly from their field- free values only for 
B* > 10^. The results in the right panels, shown as a function of B* for four values of 
u = fib/ng, show that the density dependence of this threshold value is also qualitatively 
correct. 

The upper left panel shows that there is a substantial decrease in the pressure 
associated with increasing magnetic fields for B* > 10^. This is also evident from the inset, 
which clearly shows extensive softening of the EOS. The onset of changes in the pressure as 
a function of the magnetic field may be more clearly seen in the upper right panel, in which 
results for representative densities are shown. 

The neutron effective mass m* is shown in the lower left panel, and demonstrates the 
extent to which the scalar field a is infiuenced by the presence of magnetic fields. Note that 
m* also enters in the calculation of all thermodynamic quantities. Again, it is clear that 
effects due to magnetic fields do not become significant until B* > 10^. 

The lower right panel shows that the composition of neutron-star matter changes 
significantly at high magnetic fields. The striking feature is the large increase in the proton 
fraction for B* > 10^. This has two significant effects upon the EOS. First, the protons, 
which are spin polarized, begin to dominate the contributions to thermodynamics arising 
from the baryons. This leads to a substantial softening of the EOS (see upper left panel). 
The second effect stems from the requirement of charge neutrality. Because the leptons 
provide the only source of negative charge, the lepton fraction rises commensurately with 
the proton fraction. As a result, the lepton contributions to the pressure and energy density 
are somewhat increased relative to the field-free case. However, the contributions from the 
baryons remain dominant. 

It is important to note that in order to obtain the total energy density and 
presssure relevant for neutron star structure, contributions from the electromagnetic field 
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Ef = Pf = 5V(87r) = 4.814 x 10~^5* ^ MeV fm"^ must be added to the matter energy 
density Sm and pressure P^. This has not always been done in the hterature. For B* > 10^, 
the field contributions can dominate the matter pressure, for the densities of interest, as 
shown in the upper right panel of Figure 1. 

Figure shows the dependance of H/ B on the field strength B* and the density u for 
the baseline model GM3. The so-called de Haas- van Alphen oscillations are evident and 
highlight the multi-component nature of the system. The origin of the increasing complexity 
in the oscillations may be understood by first inspecting the oscillation period when only 
a single charged species is present. When the quantity [E'j — m^)/{2qB) successively 
approaches integer values, successive Landau orbits begin to get populated resulting in 
an oscillatory structure in H/B. The width of these oscillations may be estimated by 
considering the change in magnetic field AS required to increase n-max by 1. It is found to 
be dependent upon both the strength of the magnetic field and the Fermi momenta, and is 
given by 

where -B_ and -B+ denote the fields at the beginning and end of an oscillation. In the low 
field limit 

AB^2q(^^ , (37) 

and the period goes to zero. At subnuclear densities, where muons are generally absent, 
charge neutrality forces the Fermi momenta of protons and electrons to be equal and only a 
single oscillation period exists. However, as the density increases above nuclear densities, 
the appearance of muons introduces further structure in the oscillations as a result of the 
superposition arising from each of the three charged species present. Furthermore, with 
increasing density the Fermi momenta of all particle species increase, which decreases the 
oscillation periods. The insets in each of the panels clearly show these features. 
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At large enough magnetic fields, only one Landau level is occupied, and the value of 
H/B saturates. Beyond this point, the fraction of H that the magnetization comprises 
becomes increasingly small. This is demonstrated in the lower right panel of Figure ^, in 
which the ratio H/B approaches unity for both 5* = and B* ~ 3.3 x 10^ ^ B^. However, 
for B* ~ 10^ there is a noticeable deviation of H/B from unity. Nevertheless, in all cases 
considered, \At[M./B\ does not exceed 4%, which does not represent a significant deviation 
from the case in which the magnetization is neglected. 

In Figure we compare results of the matter pressure and effective mass for the 
models HS81 (employed earlier by Chakrabarty (1996)), GMl, GM2, and ZM with a view 
towards extracting generic trends induced by strong magnetic fields. Although quantitative 
differences persist between the models, the qualitative trends of the effects of the field on 
these EOSs are shared with those of model GM3. Remaining differences are principally due 
to variations in the underlying stiffnesses, effective masses, and symmetry energies of the 
individual models. 

4. EFFECTS INCLUDING ANOMALOUS MAGNETIC MOMENTS 

We turn now to the inclusion of the anomalous magnetic moments of the nucleons. 
Johnson & Lippman (1950) first considered the inclusion of anomalous magnetic moments 
in the Dirac equation, but their formulation was noncovariant. Here, we employ the 
covariant form, suggested by Bjorken & Drell (1964), using Ci and Cn to evaluate the 
effects of magnetic fields. With the inclusion of the anomalous magnetic moments, the 
baryon field equations become 

«• [p- Qb^ + 13ml 
The derivation of the Dirac spinors is presented in Appendix A. 



^,= (E- g^^u' - gp^nj + Kk-io [l,, lA F^" • (38) 
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The energy spectrum for the protons is given by 



E 



\ 



ki + Um; 2 + 2 + i + si) q,B + sk,b\ + - ig,,/ (39) 



which may be compared with the resuh 



obtained by using the Johnson & Lippman form Kh{i/2) [y^,'y,y] F^'^ in equation ([38| ) 
for the inclusion of the magnetic moment. In both cases n and s are the principle 
quantum number and "spin" quantum number respectively. As will be shown in the 
appendix, unlike the Kb = case, the "big" components of the Dirac spinor are no longer 
eigenstates of the spin operator along the magnetic field (cr^). However, as Kp tends 
toward zero, it is clear that the proton energy spectum reduces to the expression given 



in equation (|T¥) and s corresponds to the az eigenvalue. In the nonrelativistic limit. 



k"^ <^ and 2vqpB <^ rn^ and for (upB)^ ^ m^, both of the above results reduce to 
Ep^n,s ^171* + k'^/2m* + {n + 1/2 + s/2) {qpB)/m* + SKpB, the standard nonrelativistic 
expression. 

The evaluation of the thermodynamic quantities proceeds along the lines already 
presented in §2, but with a new definition of the Fermi momentum to account for the 
presence of magnetic moments, which cause an asymmetry in the phase space in addition 
to that caused by the charged particle interactions with B. With the energy spectrum in 
equation (|39|) , 



E^- {\lm;^ + 2(n + ^ + ^s)qpB + SKpB] . (41) 



Setting 



m = \ m. 



* 2 
P 



1 1 



2[n + - + -s]qpB + SKpB , (42) 
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and substituting kf^n,s and m into the formulas for the number and energy density for the 
leptons (see § 2), one finds the analogous quantities for the protons: 

\(lp\B 



m 



(43) 
(44) 



The scalar density may be defined in terms of the energy spectrum by 

dE 



n 



dm* 



(45) 



Utilizing equation (|39D results in 
rip = J (fk 

\qp\B 

27r2 



m rrip 
m — SKpB E 



EE-;= 



m 



— SKpB 



In 



m 



(46) 



The appearance of the factor m/{m — sUpB) may be understood by inspecting the zeroth 
component of the current four-vector 



with 



7 



E 



E 



E\k=o m* + sKpB' 

where B' is taken in the rest frame of the particle. Hence, the scalar density becomes 



(47) 



(48) 



'ml + SK„B' , ^ 



d-^k- 



m rrip 
m — SKpB E 



(49) 



In a similar way, the energy spectrum of the neutrons is given by 



En,s = \lkl + (^m* ^ + kl + kl + SKnB) + ^,„^° + -QpJ 



1 







(50) 
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Including the magnetic moment for the neutron, the integral over phase space for any 
thermodynamical quantity Q may be easily evaluated by noting that at zero temperature, 
it is simply the integral over all momenta within the Fermi surface defined by 

Ef = En^s{kx, ky, kz) ■ (51) 

The integral may be written in terms of parallel and perpendicular components, 

(Q) = E ^ [ k^dk^ [ dh Q , (52) 
where a and 6 are determined by the Fermi surface to be 



E f-i V^i + m*J + SKnB , (53) 



b = y(^^^^s^«^^sj2^^7^. (54) 

With the substitution 

X = ^Jk\-\-m*^'^ + sUnB , (55) 
the integral is transformed into 

(<5) = E ( y_ y/ ^^11 <3 I - SKnSU dx J J dk\\ Qj , (56) 

where 

m = m* + SKnB . (57) 

Note that the first term is precisely the same as for the k„ = case, but with a shifted 
mass. This form for the integral over phase space is particularly useful for calculating the 
number and energy densities. Defining kf^g by 

kf,s = ^E^^-rff (58) 
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the number and energy densities take the form 



27r^ V3 



CTTl 
arcsin —— 
E] 



n 
2 



^ E 'hs + \st^nBE] 3 l^arcsin 



47r2 V2 

+ -SK„B m 

V3 " 4 



m TT 



m 



The scalar number density reads 



ml 



/c^ + m. 



* 2 I E 



which may be recast as 



ml 



/ xdx / dkii , 



Performing the integration gives 



nZ 



m 



(59) 



(60) 



(61) 



(62) 



(63) 



As in the case without magnetic moments, the pressure in beta equihbrium is given by 
P = finUb - e. 



4.1. Magnetization 



Utihzing the expressions for the energy and number densities derived above, the 
magnetization including the effects of the anomalous magnetic moments may be calculated 
using the general relation in equation ([3^). For protons and neutrons, the results are given 



by 



Mr 



B 



— r m In 



E^f + k^, 



m 



n + ^ + Is 
m — SKpB 



(64) 
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M 



n 



1 

2^ 




-SKnB m mm — > . 

.2" 3/ \ m J j 



m 3 




(65) 



The extent to which the anomalous magnetic moments alter the magnetization relative 
to the case in which they are absent may be gauged by the magnitudes of H/B in single 
component systems. For example, at fields below 2.2 x 10^^ Gauss, H/B is reduced by 
approximately 1% in a proton gas and by about 0.3% in a neutron gas. 



To assess the infiuence of the anomalous magnetic moments on the EOS, it is instructive 
to consider a charge neutral npefi gas in beta equilibrium. In addition to providing contrasts 
with the case in which only the effects of Landau quantization are considered (Lai & 
Shapiro 1991), it sets the stage for the effects to be expected for the case in which baryonic 
interactions are included. 

The magnitude of the magnetic field required to induce significant effects on the 
EOS due to the inclusion of the magnetic moments may be inferred by considering the 
field strength at which neutrons become completely polarized. From equation (|60|) , it is 



nuclear density, this leads to 5* = 1.6 x 10^. Note that this is approximately where the 
effects due to Landau quantization become large. This implies that a complete description 
of neutron-star matter in the presence of intense magnetic fields must necessarily include 
the nucleon anomalous magnetic moments. 

The equations governing the thermodynamics of the gas are simply the non-interacting 
limits of equations (|^) through ([TT]), and equation ( |5^ for the magnetization. The results 
are presented in Figure | in which the darker (lighter) shade curves show results with 



4.2. Results for the npefi Gas 



clear that complete polarization occurs when 
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(without) the inclusion of the anomalous magnetic moments. 

The left panels, in which the matter pressure is shown as functions of u and Sm, clearly 
show that the EOS is stiffened upon the inclusion of magnetic moments. For example, in 
the extreme case when the field strength approaches the proton critical field, the pressure 
is increased by an order of magnitude over the zero field case (and two orders of magnitude 
over the case in which only the effects of Landau quantization are considered) . The upper 
right panel, in which the matter pressure is shown as a function of B*, shows that above 
B* = 10^ the effects of the magnetic moments are more significant than those due to 
Landau quantization, and cannot be ignored. 

The lower right panel provides some insight into the origin of the stiffening. At field 
strengths of B* — 10^, the composition of matter is dominated by neutrons, the proton 
fraction being small, about 0.1. Neutrons, however, are spin (up) polarized due to the 
interaction of the magnetic moment with the magnetic field. With increasing S, the 
fraction of neutrons that are polarized increases leading to a corresponding increase in the 
degeneracy pressure. Upon complete polarization, this increase is halted due to the absence 
of neutrons needed to fill further spin up energy levels. This is evident from the turnover in 
the matter pressure, occuring precisely at the point when the neutrons become completely 
spin-polarized, shown in the upper right panel. 

4.3. Results for Interacting Matter 

In this section, we include the effects of baryonic interactions. Landau quantization, 
and anomalous magnetic moments. In the absence of magnetic fields, the dominant effect 
of interactions between the baryons is to substantially stiffen the EOS compared to the 
case in which interactions are omitted. This is chiefly due to the repulsive nature of the 
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baryonic interactions in beta stable matter. Notwithstanding the fact that the absolute 
magnitudes of the energy density and pressure are larger than the case in which the baryonic 
interactions are omitted, magnetic fields have many of the the qualitative effects discussed 
in the previous section. 

The results for the baseline model GM3 are shown in Figure ^, which should be 
compared with Figure |I] to assess the role of magnetic moments. The upper left panel shows 
that the stiffening of the EOS observed for the npefi gas (for B* > 10^) is also present in 
the case when interactions are included. The effects of magnetic moments are such that 
the softening caused by Landau quantiziation alone is overwhelmed, leading to an overall 
stiffening of the EOS. In fact, for fields on the order of the critical proton field, the EOS 
approaches the causal limit, Pm = ^m- As in Figure |1|, the matter pressure Pm, the effective 
mass m*, and the concentrations Yi = rii/nf, begin to differ significantly from their field-free 
values only for 5* > 10^. 

The neutron effective mass m* is shown in the lower left panel. The behavior of m* 
with B* is opposite to that shown in Figure |I|. The effects of magnetic moments cause m* 
to increase at a rate approximately equal to UpB/rrin and to become independent of density 
for B* > 10^. Note that this feature is also a consequence of complete spin polarization. 

The lower right panel shows the relative concentrations. Comparing with Figure |l], it 
is evident that the composition of matter is principally controlled by the effects of Landau 
quantization. In contrast, the stiffening of the EOS is caused primarily by terms that are 
explicitly dependent upon the magnetic moments in the pressure and energy density. 

Figure ^ to be compared with Figure ^ shows H/B as functions of both B* and u 
for the baseline model GM3. The origin of the oscillations is similar to that discussed in 
conjunction with Figure ^, but there is an overall reduction of approximately 1% in H/B 
caused chiefly by the magnetization of the neutron. 
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In Figure |^ (to be compared with Figure |^), we compare results among the models 
HS81, GMl, GM2, and ZM with the intention of extracting generic trends induced by the 
inclusion of magnetic moments. The pressure and effective masses share the qualitative 
trends exhibited by model GM3 (shown in Figure although quantitative differences 
persist between the models. The stiffness induced by the inclusion of magnetic moments 
emerges as a general trend, and remaining differences are principally due to variations in 
the underlying stiffness, effective mass, and symmetry energies of these models.. 

5. SUMMARY AND OUTLOOK 

We have developed the methodology necessary to consistently incorporate the effects 
of magnetic fields on the EOS in multicomponent, interacting matter, including a covariant 
description for the inclusion of the anomalous magentic moments of nucleons. This 
methodology is necessary because in the presence of the field all thermodynamic quantities 
inherit the dimensionful scale set by the magnetic field, which necessarily affects the 
composition and hence the EOS of matter. By employing a field theoretical-apporach which 
allows the study of models with different high density behaviors, we found that the results 
of incorporating strong magnetic fields were not very dependent upon the precise form of 
the model for the nucleon-nucleon interaction. The generic effects included softening of the 
EOS due to Landau quantization, which is, however, overwhelmed by stiffening due to the 
incorporation of the anomalous magnetic moments of the nucleons. These effects become 
significant for fields in excess of B* ~ 10^, for which neutrons become completely spin 
polarized. Note that this field strength is substantially less than the proton critical field. In 
addition, the inclusion of ultra-strong magnetic fields leads to a reduction in the electron 
chemical potential and an increase in proton fraction. These compositional changes have 
implications for neutrino emission via the direct Urea process and, thus, for the cooling 
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of neutron stars. The magnetization of the matter never appears to become very large, 
as the value of \H/B\ never deviates from unity by more than a few percent. However, it 
remains to be seen what effects the magnetization of matter will have on the structure and 
transport properties of neutron stars. 

It is worthwhile to note here that the qualitative effects of strong magnetic fields 
found in the relativistic field-theoretical description of dense matter would also be found in 
non-relativistic potential models. This is because the phase space of charged particles is 
similarly affected in both approaches by the presence of magnetic fields. The effects due 
to the anomalous magnetic moments would, however, enter linearly in a non-relativistic 
approach (see §4), and would thus be more dramatic in this case. It would be also be 
instructive to study the effects of magnetic fields including many-body correlations. 

It would be useful to also consider cases in which strangeness-bearing hyperons, a 
Bose (pion or kaon) condensate or quarks, are present in dense matter. The covariant 
description of the anomalous magnetic moments developed in this work may be utilized to 
include hyperons, which are likely to be present in dense matter (Glendenning 1982, 1985; 
Weber & Weigel 1985; Kapusta & Olive 1990; EUis, Kapusta & Olive 1991; Glendenning 
& Moszkowski 1991; Sumiyoshi & Toki 1994; Prakash et al. 1997 and references therein). 
The anomalous magnetic moments of hyperons are mostly known. The negatively charged 
hyperons, the neutral A, and S° all have negative anomalous magnetic moments. and TP 
are the only hyperons with positive anomalous magnetic moments. The effects of Landau 
quantization on hyperons would be to soften the EOS relative to the case in which magnetic 
fields are absent. However, in the presence of strong magnetic fields, all of the hyperons 
will be spin polarized due to magnetic moment interactions with the field. This would 
cause their degeneracy pressures to increase compared to the field-free case. The resultant 
of these two opposing effects will depend on the relative concentrations of the various 
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hyperons, which in turn depends sensitively on the hyperon-meson interactions for which 
only a modest amount of guidance is available (Glendenning & Moszkowski 1991, Knorren, 
Prakash & Ellis 1995, Schaffner & Mishustin 1996). For choices of E~— meson interactions 
that favor the appearence of E~ hyperons at relatively low densities, the concentrations of 
the positively charged particles, p and E+, may be expected to increase in the presence of 
strong magnetic fields. It would thus appear that the effects of including hyperons will not 
drastically alter the qualitative trends of increasing the concentrations of positively charged 
particles found in the case of npe^ matter. The main physical effects found in the absence 
of hyperons, namely increasing the stiffness of matter, and allowing the direct Urea process 
(Lattimer et al. 1991; Prakash et al. 1992) to occur, probably would not change, either. 
Feedback effects due to mass and energy shifts may, however, alter these expectations. 
Thus, detailed calculations are required to ascertain the influence of magnetic fields in 
multi-component matter. Work on this topic is currently in progress and will be reported 
separately. 

It is intriguing that Bosons (pions and kaons), which have zero magnetic moment, 
do not feel the magnetic fields as fermions do. Similarly, quarks without sub-structure 
also have no anomalous magnetic moments. Thus, intense magnetic fields in the cores of 
stars containing a Bose condensate or quark matter might serve as a useful discriminant 
compared to those containing baryonic matter. 

Work is in progress (Cardall et al. 1999) to complete a fully self-consistent calculation 
of neutron star structure including the combined effects of the direct effects of magnetic 
fields on the EOS, which we have developed in this paper, and general relativistic structure. 
The findings will help answer questions concerning the largest frozen-in magnetic field that 
a stationary neutron star can possess, and what the structure of stars with ultra-strong 
fields might be. It must be borne in mind, however, that for super-strong fields (much 
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higher than B^, which is the highest field considered in this work), the energy density in 
the field would be significantly higher than the baryon mass energy density. Under such 
conditions, the internal structure of the baryons will be affected and alternative descriptions 
for the EOS will become necessary. 

We thank Hans Hansson for constructive suggestions concerning the covariant 
description of the anomalous magnetic moments. This work was supported in part by the 
NASA ATP Grant # NAG 52863, and by the USDOE grants DOE/DE-FG02-87ER-40317 
& DOE/DE-FG02-88ER-40388. 
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A. SPINORS AND ENERGY SPECTRA FOR BARYONS WITH 
ANOMALOUS MAGNETIC MOMENTS 

In this appendix, we derive relations for the spinors and energy spectra for baryons 
with anomalous magnetic moments. The Dirac equation is 



a* 



(p- q^A) + 13ml + pa.KbB] = ^o,f.^fe , (Al) 



where the effective momentum is given by vf = p — qbA and Hb denotes the baryon anomalous 
magnetic moment. The energy Eq i, denotes the baryon energy eigenvalues when the meson 
fields are absent and are related to the neutron and proton energy spectra given in equations 
(0) and (|D|) by 

En,s = Eo,n + 9u.,UJ^ + ^9p,P^ (A2) 

Ep,n,s = Eo^p + g^^uj^ - ^gp,p° , (A3) 

repectively. Separating \l'b in to "big" and "small" components, we obtain 

{Eo^b - ml - KbBaz) = (a'vr) x (A4) 
{Eo^b + ml + KbBa,)x = (a-7r)0. (A5) 

Writing x in terms of (taking care to note that the terms on the left hand side of these 
equations are no longer proportional to the identity matrix because of the presence of the 
magnetic moments), equation (Xi) becomes 



(Eo,fe - - KbBa^) (f) = (a-vr) — ■ — {a-ii) <p . (A6) 

(Eo,6 + ml) - (KbB) 



Note that the term with cxz does not commute with the momentum operators. Therefore, 

[Eo^b + ml) - {KbB) [Eo^b + ml) - {KbB) 

(A7) 



- 30 - 



This may be rewritten as 

Fs4> = {S'T^f (p-as (a«7f) 7r^0 , (A8) 
where Fg and are defined as 

Fs = {Eo^b - K.bBa^f - ml ^ 

(3o — 7^ 7^ — . (Ay) 

iEo,b + mlf-{KbBf ^ ' 

At this point it is necessary to consider individually the cases of the protons and neutrons. 



Protons 



qpB 



TTy. 



The fact that [7ra;,7ry] = ih{qB/c) suggests the transformations 

Pi -- 

Using the identities 

(a»a) {(f'bj = a*b + ia* (a xbj , vf x vf = iqpBz 



and the above transformations, equation (|A8D becomes 
Fs(f) = [qpB (pi + - a^) + P^] - a. J ^ 



(AlO) 



(All) 



(A12) 



The similarities between equation (|A12|) and that for the leptons (see, for example, Itzykson 
& Zuber 1984), suggests the ansatz for the spin up spinor 



0+1 = e*''"^"~ 



/ 



(A13) 



\ 'i'^p,n,+lHn+l{i) 

The two coupled differential equations for the components of (equations ( |A12D ) reduce to 
two coupled algebraic equations for the eigenvalues of i?o,p and ujp^n,+i- Explicitly, 



F+i = 2{n+l)qpB + {I- a+i)k^- a+i2{n + l)JqpBk^ujp^n,+i 



F_i = 2{n + l)qpB + {1 + a-i)kl- a^ik^JqpBujJ^ j^^. 



(A14) 
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These may be solved to give 



Eo,p,+i = y A;2 + (^^ml'' + 2{n + l)qpB + n^B^ ' 



(A15) 



With this result, it is straightforward to solve for a;p,„,+i. Lacking a simple expression, we 
shall continue to refer to it as u}p^n,+i- Substituting this solution for 0+i into the expression 
for X gives the Dirac spinor 



K^i = Ne 



Hn-l{i) 



\ 



Eo,p,+i+m*-\-KpB 
Eo,p,+i+'m;-KpB -"nUj / 



A similar method may be employed to find an ansatz for the spin down spinor. 



ik z ^ i^P,n,-iHn-i{i) 
0_i = e 2 



with the energy eigenvalue 
and the Dirac spinor 



ml 2 + 2nqpB - k^B^ , 



V 



-2/1/ ^ <ii,B-i^i,,n.-ik 



Eo,p,-i + m;+KpB ^n-liO 
Eo,p,-i+m*p-KpB-'^n\^) I 



(A16) 



(A17) 



(A18) 



(A19) 



While all quantities in this work have been calculated in the zero temperature 
approximation, requiring only the postive energy spinors, for completeness the negative 
energy Dirac spinors are presented below. For the protons these may be determined in 
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much the same manner as that employed for the positive energy spinors. Defining 



\2 * 2 



2KpB (Eq^p - m* + KpBa, 



(A20) 

the equation for x takes the same form as equation ([A8|) where F~ and a~ replace Fs 
and respectively. As a result, precisely the same formalism employed to determine the 
positive energy spinors may be used to determine the negative energy spinors. The Dirac 
spinor corresponding to the energy eigenvalue 



0,p,+l 



kl + Um* 2 + 2 + 1) gp5 + KpB 



(A22) 



ra,+l 



V 



(A23) 



is given by 



where is defined by replacing Fs and in equations ( |A14] ). Similarly, the Dirac 

spinor corresponding to the energy eigenvalue 



0,p,-l 



A;2 + Jm* 2 + 2nqpB - KpB 



(A24) 



is given by 



/ -2inU qpB-iu) ^ _-^kz 



^ n, — 1 



A^e" 



p,n, — 1 



p,n-lHn-l{C) 



(A25) 
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Neutrons 

In this case, the trial wave function has the same form as the free particle solutions 
with unknown coefficients, which may be determined in a manner analougous to that 
employed for the protons. Define 

G^ = Fs-k'^ + a,a,kl . (A26) 

Then, equation (|A8|) becomes 

Gs(j) = -as{a^k^ + ayky)k^(j) . (A27) 

Note that both Gg and are diagonal and therefore the off-diagonal terms have been 
isolated on the right-hand side of equation ( [A27| ). The similarities with the case in which 
Kn = 0, namely the quadratic nature of the momentum operators, suggests the form 

= e-'*='"- ( ) ' ^"^^^^ 
Using equation ( |A28| ) , we obtain the coupled algebraic equations 

G^iu = —a^iikx — iky)kzV 

G-iv = —a-i{kx + iky)kzU. (A29) 

Combining these gives 

G+iG.i = a+ia_i{kl + kl)kl , (A30) 
which may be solved for the energy eigenvalue 



E,,n,s = \jkl + (^m; 2 + ^2 + ^2 + ,^^B) . (A31) 



The eigenvectors may be determined, up to a normalization, by setting 



u = 1 ^ V = — {kx + iky)kz (A32) 



u 



a_i 

g7. 



{k,+lky)k, (A33) 
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in equation ( |A29|) . It is clear from direct substitution that the first gives the s = +1 and 
the second the s = — 1 spinors. Inserting these into equation ( [A28| ) and then into the 
expression for x gives the neutron Dirac spinors 



= Ne 



1 



G. 



{kx + iky) 



Eo,n, + l+M*+K„B 



\ Eo,n, + l+M*-K„B ' / 



(A34) 



1+1 



+1 



G. 



1 G+i^^ 



iky) kz 



Eo,n,-l+M*+K.„B 



V Eo,„-i+M*-K„B / 



(A35) 



In order to determine the negative energy Dirac spinors for the neutron, an approach 
analogous to that employed in determining the negative energy Dirac spinors for the protons 
may be used. Define Gj by replacing Fs and by and a~ , respectively, in equation 
( |A26| ). As in the case of the protons, this produces an equation for x which is of the same 
form as that employed for in the derivation of the positive energy spinors. Proceeding in 
the same manner as before, one finds that the Dirac spinors corresponding to the energy 
eigenvalues 

(A36) 



0,n,s 



kl + + kl + kl + SKnB 



-35- 



are given by 



= Ne 



-1 ^ 



1- 



Eo,n,+l-m*+KnB 
1 



(A37) 



Ne 



- ~l~ iky) kz 



l-^kl 

G7, 



Eo^n.-l-m^~HinB 
+1 , 



Eo,n,-i-m* +re„S 



+1 



(ASS) 
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TABLE 1 

NUCLEON-MESON COUPLING CONSTANTS 



Model 


ris 


-B/A 


M*/M 




^sym 








h 


c 


HS81 


0.148 


15.75 


0.54 


545 


35.0 


3.974 


3.477 


2.069 


0.0 


0.0 


GMl 


0.153 


16.30 


0.70 


300 


32.5 


3.434 


2.674 


2.100 


0.002947 


-0.001070 


GM2 


0.153 


16.30 


0.78 


300 


32.5 


3.025 


2.195 


2.189 


0.003478 


0.01328 


GM3 


0.153 


16.30 


0.78 


240 


32.5 


3.151 


2.195 


2.189 


0.008659 


-0.002421 


ZM 


0.160 


16.00 


0.86 


225 


32.5 


2.736 


1.617 


2.185 


0.0 


0.0 



NOTE.- Coupling constants for the HS (Horowitz & Serot 1981), GMl-3 (Glendenning & 
Moszkowski 1991), and ZM (Zimanyi & Moszkowski 1990; 1992) models. The couphngs 
arc chosen to reproduce the binding energy B/A (MeV), the nuclear saturation density 
(fm~^), the Dirac effective mass M* in units of the baryon mass M, and the the symmetry 
energy a^y^{MeV). The nuclear matter compression modulus Kq (MeV) for the different 
models are also listed. 
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FIGURE CAPTIONS 

FIG. 1.- Matter pressure P^, nucleon Dirac effective mass m*^/mn, and concentrations 

= ni/rii, as functions of the density u = nb/us (left panels; Ug = 0.16 fm~^ is the 
fiducial nuclear saturation density) and magnetic field strength B* = B / (right panels; 

= 4.414 X 10^^ Gauss is the electron critical field), for the model GM3. The inset in the 
upper left panel shows Pm as a function of the matter energy density Sm- The curve labeled 
Pf in the upper right panel shows the B'^/Stt contribution to the total pressure. The inset 
in the lower left panel shows the effective mass as a function of B*. In the lower right panel, 
the electron and neutron concentrations have been suppressed for clarity {Ye — Yp — Yfj^ and 

Yr,^l-Yp). 

FIG. 2.- The ratio of the induced to apphed magnetic field H/B as functions of the density 
and magnetic field strength, for the model GM3. The insets show H/B in expanded scales 
to highlight the effects of including several components. 

FIG. 3.- Matter pressure Pm and the nucleon Dirac effective mass m*^/mn for the models 
shown in Table 1 (with the exception of model GM3, whose results are displayed in Figures 
1 and 2), as functions of the density and magnetic field strength. The insets in the left 
panels show Pm as a function of the matter energy density Sm- 

FIG. 4.- Matter pressure Pm and concentrations Yi = rii/rib as functions of the density and 
magnetic field strength for a charge neutral, beta-equilibrated, non-interacting npe/j, gas 
with and without the inclusion of the nucleon anomalous magnetic moments Kb. The curve 
labeled Pf in the upper right panel shows the B'^/Sn contribution to the total pressure. 
The lower left panel shows the enhancement in the pressure, as a function of energy density, 
due to the presence of magnetic fields. In the lower right panel, the electron and neutron 
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concentrations have been suppressed for clarity {Yg — Yp — Y^ and Y^ — 1 — Y^ . 

FIG. 5.- Same as Figure 1, except that the nucleon anomalous magnetic moments are now 
included. 

FIG. 6.- Same as Figure 2, except that the nucleon anomalous magnetic moments are now 
included. 

FIG. 7.- Same as Figure 3, except that the nucleon anomalous magnetic moments are now 
included. 




Fig. 1.^ 
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Fig. 2.- 




Fig. 3.- 
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Fig. 6.— 




Fig. 7.- 



